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ABSTRACT 


A mathematical  model  of  a zonally  uniform  ocean  is 
constructed  in  order  to  investigate  the  steady  circulation  and 
the  nonhomogeneous  distribution  of  density  maintained  by 
elementary  thermodynamic  processes.  The  thermodynamic  processes 
are  considered  implicitly  by  introducing  an  equation  to  represent 
the  transfer  of  mass  by  convection  and  diffusion.  The  thermo- 
dynamic processes  which  affect  the  density  in  the  surface  layers 
of  the  ocean  are  replaced  by  a distribution  of  equivalent  mass 
flux  across  the  surface, 

A simplified  system  of  equations  is  derived  from  the 
general  system  of  hydrodynamic  equations  representing  the  model 
by  neglecting  small  terms  in  the  equations.  The  simplified 
system  of  equations  is  nonlinear  and  of  the  boundary-layer  type. 
An  approximate  analysis  is  carried  out  using  simplified 
boundary- layer  methods  and  the  following  types  of  flow  with 
their  associated  density  fields  are  studied  in  detail:  (A)  the 
flow  set  up  in  a ' omogeneous  ocean  by  surface  wind  stress,  (B) 
the  flow  and  the  density  field  set  up  by  a given  mass  flux 
across  the  ocean  surface  in  the  absence  of  wind  stress,  (C)  the 
flow  and  the  density  field  set  up  by  a given  mass  flux  across 
the  ocean  surface  in  the  presence  of  wind  stress.  For  purposes 
of  comparison,  the  same  distributions  of  wind  stress  and  surface 
mass  flux  are  used  in  studying  the  three  different  types  of 
flow. 

The  different  types  of  flow  are  Illustrated  by  diagrams 
showing  the  streamlines  of  the  meridional  circulation  and  the 
distribution  of  isopycnals  in  a meridional  section.  The  vertical 
distributions  of  the  zonal  velocity  component  for  the  different 
types  of  flow  are  also  compared. 

From  the  analysis  of  this  model,  it  is  concluded  that 
thermodynamic  processes  affecting  the  density  of  sea  water  cannot 
be  neglected  in  the  study  of  zonally  uniform  currents  in  the 
ocean.  Thus,  it  does  not  appear  to  be  possible  to  account  for 
the  observed  features  of  an  ocean  current,  such  as  the  Antarctic 
Circumpolar  Current,  in  terms  of  surface  wind  stress  alone. 
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1 . Introduction 

The  dynamics  of  atmospheric  and  oceanic  circulation  has 
been  studied  extensively  in  the  past  with  the  result  that  many 
theoretical  models  based  on  the  hydrodynamic  equations  have  been 
constructed  to  interpret  various  aspects  of  the  balance  of  forces 
achieved  in  large-scale  fluid  motion.  The  3tudy  of  tho  thermo- 
dynamic processes  by  which  some  of  these  forces  are  created  and 
maintained  has  not  been  very  successful  from  a theoretical  point 
of  view  because  of  the  difficulties  involved  in  treating  these 
processes  analytically.  As  a consequence,  discussions  of  the 
effects  of  thermodynamic  processes  on  large-scale  fluid  motion 
have  been  qualitative  in  nature.  However,  an  adequate  under- 
standing of  atmospheric  and  oceanic  circulation  cannot  be  achieved 
until  tho  baroclinic  structure  of  the  atmosphere  and  ocean  can  be 
interpreted  quantitatively  in  terms  of  suitable  models  which 
include  thermodynamic  processes. 

The  purpose  of  the  present  investigation  is  to  consider 
a theoretical  model  of  a zonally  uniform  ocoan,  i.e.,  an  ocean 
having  properties  that  do  not  vary  with  longitude.  Thermodynamic 
procossos  are  introduced  in  a sufficiently  simple  manner  to  allow 
analytical  treatment  of  the  model. 

Many  Investigations  of  atmospheric  and  oceanic  circulation 
have  been  made  using  models  based  on  zonally  uniform  flow,  so  that 
a considerable  amount  of  literature  on  the  dynamics  of  zonally 
uniform  flow  5s  available.  The  classical  Kadley  model  of  atmos- 
pheric circulation  is  an  example.  It  is  related  to  the  present 
model  In  that  the  qualitative  effects  of  thermodynamic  processes 
are  considered,  Rossby  ( 1 91+1  ) has  described  the  basic  concepts 
as  well  as  the  limitations  of  the  Hadley  model.  Attempts  have 
been  made  to  construct  a quantitative  Hadley  model.  Dorodnitsyn, 
Izvekov,  and  Schwetz  (1939),  for  example,  constructed  a semi- 
empirical  model  of  the  Hadley  type  using  observed  distributions 
of  temperature  and  surface  pressure  in  the  atmosphere;  however, 
they  avoided  explicit  discussion  of  the  thermodynamic  processes 
which  accompany  the  resulting  atmospheric  circulation. 
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The  thermodynamic  processes  occurring  in  the  ocean  are 
recognized  to  be  of  importance  in  their  effect  on  the  circulation, 
particularly  in  the  polar  regions.  Deacon  (1933)  and  Sverdrup 
(1933)  indicate  that  thermodynamic  processes  play  an  important 
role  in  determining  the  behavior  of  the  Antarctic  Circumpolar 
Current  (ACC ) . 

The  ACC  is  tho  prototype  of  the  present  model  in  that  the 
dimensions  and  the  range  of  variation  of  the  various  quantities 
considered  in  the  model  are  chosen  to  be  of  the  order  of  magnitude 
observed  in  the  ACC.  However,  the  ACC  is  not  zonally  uniform.  It 
has  been  suggested  by  Hunk  and  Palmen  (1951)  that  the  departure 
from  zonal  uniformity  may  be  necessary  for  a dynamic  equilibrium 
to  be  achieved.  As  all  zonal  variations  are  excluded  in  the 
present  model,  the  results  obtained  probably  have  limited  applica- 
bility to  actual  conditions  existing  in  the  ACC. 

2.  General  problem 

The  problem  considered  here  is,  in  its  most  general  sense, 
that  of  determining  the  circulation  and  tho  distribution  of  physi- 
cal and  chemical  properties  such  as  density,  temperature,  salinity, 
etc.,  produced  by  a given  surface  wind  stress  and  a given  flux  of 
heat  and  water  across  the  ocean  surface.  Only  a small  portion  of 
the  general  problem  Is  actually  considered  and  many  simplifying 
assumptions  must  be  made  before  explicit  results  can  be  obtained. 

At  least  two  fundamental  difficulties  are  encountered  in 
setting  up  a theoretical  model  to  include  thermodynamic  processes, 
and  these  difficulties  are  avoided  rather  than  resolved  in  the 
present  investigation.  The  first  difficulty  is  encountered  in 
setting  up  the  dynamic  portion  of  the  model.  Theoretically,  one 
might  expect  that  the  torque  of  the  zonal  component  of  surface 
wind  stress  is  either  transmitted  vertically  to  the  ocean  bottom 
or  laterally  to  the  continental  coasts.  If  the  torque  is  trans- 
mitted to  the  ocean  bottom,  there  are  strong  currents  near  the 
bottom,.  There  is  no  conclusive  observational  evidence  to  indicate 
whether  or  not  such  bottom  currents  exist.  Munk  and  Palmen  (1951) 


say  that  lateral  stresses  in  the  ACC  do  not  appear  to  be  large 
enough  to  balance  the  wind  stress.  They  suggest  that  tho  balance 
could  be  accomplished  by  the  "mountain  offoct"  of  submarine  ridges 
which,  in  effect,  produce  a retarding  torque  on  the  ACC.  The 
"mountain  effect"  cannot  be  incorporated  into  the  present  model 
except  by  postulating  ad  hoc  a zonal  force  of  an  arbitrary  nature. 
The  present  model  is  set  up  so  that  the  zonal  component  of  wind 
stress  is  transmitted  almost  entirely  to  the  ocean  bottom* 

The  second  difficulty  is  encountered  in  introducing 
thermodynamic  processes.  The  important  thermodynamic  processes 
which  affect  the  density  of  the  water  act  in  the  surface  layers 
of  the  ocean.  It  has  not  been  found  possible  to  consider  processes 
which  tend  to  increase  the  density  near  the  surface  of  the  ocean. 

If  such  processes  are  considered,  the  resulting  density  distribu- 
tion is  unstable,  i.o.,  the  density  increases  upward,  in  the  region 
where  these  processes  occur.  This  difficulty  is  avoided  by  consid- 
ering only  those  thermodynamic  processes  whose  net  effect  is  to 
lower  the  density  near  the  ocean  surface  and  by  limiting  tho  dis- 
cussion to  an  interior  region  of  the  ocean  away  from  continental 
coasts.  It  Is  possible  then  to  consider  a steady  state  In  a lim- 
ited region  of  the  ocean  by  regarding  the  ocean  outside  of  this 
region  as  a reservoir  which  can  supply  water,  heat,  salt,  etc., 
at  a rate  necessary  to  maintain  the  steady  state. 

3 . Simplifying  assumptions 

The  basic  assumptions  necessary  to  set  up  the  mathematical 
equations  of  tho  model  are  given  below.  More  specific  assumptions 
are  discussed  as  they  arise  in  the  analysis  of  the  model. 

The  ocean  basin  is  taken  to  be  zonally  uniform  and  of  con- 
stant depth  with  a level  bottom.  Lateral  boundaries  are  not  ex- 
plicitly considered.  The  coordinates,  of  a point  in  the  ocean  can 
be  specified  by  the  latitude,  longitude,  and  height  above  the 
ocean  bottom  measured  along  the  local  vertical.  It  is  assumed 
that  this  coordinate  system  can  be  replaced  with  sufficient  accu- 
racy by  a spherical  coordinate  system  (X,  $ , r)  in  which  gravity 
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acts  in  the  radial  direction.  Hero,  K is  the  longitude,  / tho 
latitude,  and  r tho  distance  of  the  point  from  the  center  of  the 
earth.  The  distance  r is  approximated  by  a + z,  whore  a is  the 
mean  radius  of  the  earth  and  z is  the  height  above  the  ocoan 
bottom. 

The  velocity  components  relative  to  the  earth  (u,  v,  w) 
aro  chosen  so  that  u is  the  eastward  component,  v the  northward 
component,  and  w the  upward  component. 

Because  the  density  variations  in  the  ocean  aro  a small 
fraction  of  the  total  density,  tho  variations  are  neglected  in 
the  equations  of  motion  except  in  the  torm  representing  tho  force 
of  gravity,  Tho  compressibility  cf  the  water  is  neglected  every- 
where, so  that  a stream  function  can  bo  introduced  for  the  circu- 
lation in  the  meridional  plane.  The  flux  of  wator  across  the 
ocean  surfaco  is  negligiblo  as  far  as  the  stream  function  is 
concerned,  because  the  amount  of  water  exchanged  with  the  atmos- 
phero  by  evaporation  and  precipitation  is  small  in  comparison  with 
the  amount  of  water  transported  by  the  meridional  currents. 
Howovor,  the  effect  of  tho  flux  of  water  on  the  density  cannot  be 
neglected  because  the  water  added  is  generally  of  a different 
density.  In  its  effect  on  the  density,  the  flux  of  water  can  be 
interpreted  as  a flux  of  mass  numerically  equal  to  the  difference 
in  mass  between  tho  volume  of  water  added  across  a unit  aroa  of 
ocean  surface  per  unit  time  and  an  equal  volumo  of  oceanic  water. 
Similarly,  the  volumo  change  due  to  thermal  expansion  is  not 
significant  in  comparison  with  the  total  volume  of  wator  in 
circulation  but  the  density  change  due  to  the  thermal  expansion 
is  dynamically  significant.  The  heat  flux  must  be  included  in 
terms  of  an  equivalent  mass  flux  across  the  surface  of  the  ocean. 

Tho  diffusion  and  convection  of  heat,  salt,  etc.,  in  the 
interior  of  tho  ouean  aro  not  takon  into  account  separately,  but 
their  effect  is  taken  into  account  by  introducing  an  equation  to 
represent  the  diffusion  and  convection  of  mass.  The  assumptions 
under  which  this  simplification  is  justified  are  given  in  Appendix 
I,  which  includes  the  derivation  of  the  mass- transfer  equation. 

Instead  of  using  nonisotropic  eddy  coefficients  having 


- 5 - 


large  values  to  characterize  horizontal  transfers  of  momentum  and 
mass,  the  nonlinear  terms  representing  convective  transfers  of 
these  quantities  are  retained  in  the  equations  and  the  diffusive 
processes  are  assumed  to  be  characterized  by  constant,  isotropic 
kinematic  eddy  coefficients.  The  kinematic  eddy  coefficients  of 
viscosity  and  diffusivity  are  assumed  to  be  of  different  orders 
of  magnitude  for  reasons  given  in  Appendix  II. 

A basic  assumption  underlying  the  entire  development  of 
the  model  is  that  the  circulation  produced  by  steady,  zonally 
uniform  distributions  of  wind  stress  and  mass  flux  is  steady  and 
zonally  uniform.  The  dynamic  stability  of  the  circulation  is, 
therefore,  not  investigated  and  the  flow  is  treated  as  being 
stable « 

I4..  Equations  of  motion 

On  tho  basis  of  the  preceding  assumptions,  tho  equations 
governing  the  steady,  zonally  uniform  motion  and  tho  distribution 
of  density  in  the  ocean  may  be  writton 
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is  i a 2 

(cos  / v)  + — 2 — (r  w)  = 0 

r cos  <&  3^  r cr 


v dp  ap  p 

+ w — = K v p,  ( 

r ctf  cr  H 

where  fi  is  the  earth’s  angular  velocity,  p the  density  considered 
as  a function  of  the  coordinates,  p a constant  mean  density,  and 
p the  pressure.  The  kinematic  eddy  coefficients  of  viscosity  and 
diffusivity  are  designated  by  K and  Kp  respectively,  and  gravity 
by  g.  The  Laplacian  in  this  coordinate  system  is 

2 1 3 / £ \ 1 c. ' £ s \ 

<7  = — *j COS  / — ',+  ~~2  — jr  — j. 

r cos  yi  ddl  r dr  V dr/ 

A stream  function  t is  introduced  for  the  meridional 


circulation  so  that 
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The  quantities  M,  x,  D , and  G are  introduced  and  defined  as 
follows : 

M = r cos  ft  u, 

the  unit-mass  angular  momentum  of  the  zonal  flow, 

x = sin  ft 
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Tho  differential  operators  D and  G arc  commutative.  The  operator 
G represents  differentiation  parallel  to  the  earth's  axis  of 
rotation. 

Equations  (1)  to  (5)  can  be  expressed  more  compactly  in 
torms  of  those  quantities  as  follows: 
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is  the  Jacobian  of  tho  functions  ^ and  f with  respect  to  r and  x. 
Elimination  of  the  pressure  from  (7)  and  (8)  yields 
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The  boundary  conditions  to  bo  satisfied  along  the  bottom 
of  the  ocean,  r = a,  are  that  all  components  of  the  velocity  are 
zero  and  that  the  mass  flux  through  the  bottom  is  zero.  Along 
the  surfaco  of  the  ocean,  r = a + h,  the  boundary  conditions  are 
that  the  surface  stress  must  equal  the  ^iven  wind  stress  and  that 
the  mass  flux  across  the  surfaco  must  equal  the  given  mass  flux. 

Further  restrictions  must  be  made  in  order  for  a steady 
state  to  exist.  The  steady=-s tato  restrictions  are  that  there  can 
be  no  net  volume  transport  across  any  latitude  circle  and  that 
the  net  flux  of  angular  momentum  and  of  mass  into  the  ocean  must 
be  zero. 

The  boundary  conditions  may  bo  expressed  mathematically 
as  follows: 

along  the  ocean  bottom,  r = a, 

at/3*  = = 0 

M = 0 


dp/dr  - 0, 


and  along  the  ocean  surfaco,  r = a + h, 

aijr/ax  = o 


P a m 

Kr n 

8r  r 


3p  1 


Here,  n:  and  % are  the  meridional  and  zonal  component  respectively 
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of  the  surface  wind  stress,  T is  the  torque,  about  the  earth’s 
axis  of  rotation,  produced  by  tho  surface  wind  stress,  and  Q is 
the  mass  flux  across  tho  surface  of  the  ocoan. 

Tho  steady-state  restrictions  on  T*,  and  Q*  are 

jftnda  = jJT*nda  = J|q*nda  = 0, 

where  T and  Q,  arc  the  zonal  torque  and  mass  flux  respectively s 
The  surface  integrals  are  taken  over  the  entire  ocean  boundary, 
with  n denoting  a unit  vector,  normal  to  the  boundary,  directed 
into  the  ocean.  As  the  boundary  of  the  ocean  must  be  a stream- 
line, the  integral  condition  on  ljl  may  be  replaced  by  the  condi- 
tion ijj  = 0 along  the  boundary.  If  the  mass  flux  through  the 

•ft 

ocean  bottom  is  zero,  the  condition  on  q , for  an  ocean  lying 
between  x^  and  x^,  reduces  to 

2qdx  =••  o . 

Thus,  if  the  surface  mass  flux  is  different  from  zero,  it  must 
have  both  positive  and  negative  values  to  satisfy  the  steady- 
state  condition. 

5.  Dimensional  analysis 

The  mathematical  problem  set  up  in  the  previous  section 
is  too  complex  to  solve,  so  that  further  simplification  is  neces- 
sary, This  simplification  is  achieved  by  neglecting  terms  which 
are  estimated  to  be  small  in  the  equations  representing  the  model. 
An  estimate  of  the  relative  size  of  the  various  terms  is  made  by 
assigning  characteristic  magnitudes  to  the  dimensional  quantities 
entering  into  the  equations . 

Accurate  estimates  can  be  made  of  the  earth’s  mean  radius 
a,  the  earth’s  rate  of  rotation  si,  the  mean  depth  of  the  ocean  h, 
gravity  g,  and  the  mean  density  ]d.  The  magnitudes  of  the  surface 
wind  stress  t,  the  surface  mass  flux  Q,,  and  the  fractional  range 
of  density  variation  y can  be  estimated  roughly  from  observations 
taken  in  the  ocean,  ho  satisfactory  estimates  are  available  from 
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observations  as  to  the  magnitudes  of  the  kinematic  eddy  coeffi- 
cients K and  Kh. 

P 

The  equations  representing  the  present  model  can  best  be 
studied  by  transforming  them  into  nondimens ional  form.  By 
studying  the  nondimens ional  form  of  the  equations,  it  is  possible 
to  evaluate  the  effects  of  both  large  and  small  values  of  the 
eddy  coefficients. 

The  dependent  and  independent  variables  are  converted 
into  nondimens ional  form,  denoted  by  primes,  through  the 
following  transformation: 

r = aq'  = a(1  + az»),  a = h/a 

x = x 1 

p = p(1  + Yp') 

Tjr  = • 

M = aUM » 

T = aT^T » 


Q = 

where  •'4r  and  u are  characteristic  values  of  the  meridional  stream 
function  and  the  zonal  velocity  component  respectively,  and 
are  maximum  values  of  the  wind  stress  and  surface  mass  flux 
respectively. 

The  equations  obtained  after  performing  the  transfor- 
mation have  the  form 
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where 


1 3(t«,p') 

a2h  q»2  8 (z  * ,x ' ) 


K 


- -%  V 


» rs  » 


a 1 - 

G 1 ~ x * + a 

az»  q» 


2 

x'^  a 


ax» 


(13) 


and 


D»‘ 


. .2  -.2 

1 - x1  a 


,2  = 


3z 


1 


+ a 


— -£ 

q»^  3x«c 


-,) + a' 


(1  - xc) 


q 1 f 3z  ’ \ 3z  » / 3x » 

The  boundary  conditions,  in  nondimens ional  form,  are 

0 ~ 

0 


A 

axj  \ 


= at’/az- 
M* 

3p»/3z> 


= 0 


J-L'*±)  = 

az « \ az »/ 


> at  z * = 0, 


t ’ = 0 


ah2T 


^ q.(1  - x.2)*  ,i 

PK-i' 

2 i 

tt’ 

m 


L at  z » = 1 


= (>y£)q'(i  - x»  )2  t 


» t 

s 


L 


[ 

\ 
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2 ^ ( - 2 

(q'  M» 

3z»  V / 


= XT'  ;>  at  z ' = 1 . 


8p ' 
cz » 

Here, X,  , and  £ represent  nondimens ional  constants. 

The  nonlinear  term  in  (11)  is  of  significance  only  if  the 
zonal  velocity  component  is  a considerable  fraction  of  the  equa- 
torial velocity Xla  of  the  earth's  surface.  As  such  high  veloci- 
ties are  not  present  in  the  ^cean,  equation  (11)  must  represent  a 
balance  essentially  between  the  Coriolis  force  and  the  gradient 
of  the  zonal  stress  component.  In  order  to  express  this  balance 
by  nondimens ional  terms  of  equal  magnitude,  the  dimensional  coef- 
ficients of  (11)  are  chosen  so  that 


h^H 

ypK. 


= 


2ify/h  = KaU/h2. 

Similarly,  it  can  be  shown  that  (12)  represents  a balance  mainly 
between  the  horizontal  gradient  of  the  force  due  to  gravity  and 
the  vertical  gradient  of  the  Coriolis  force,  i.e.,  the  zonal 
current  is  approximately  gcostrophic.  The  balance  is  expressed 
in  nondimens ional  form  by  choosing 

2fiaU/h  = Yg  • 

The  characteristic  values  of  \|/  and  U are,  therefore, 

■’4r  « YgK/  ( 2i2) 2 


U = Yga /2£l. 

The  characteristic  value  of  y is  chosen  so  that  = 1,  i.e., 

Y = hQM/pKp. 

Equations  (10),  (11),  and  (12)  can  be  written  in  nondimen- 


I 
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sional  fom  using  the  nondimens ional  parameters  6,  and  ja 
where 


5 = &J/na 


€ = K/2fthc 


= 6K/Kp. 

The  nonaimens ional  equations,  with  primes  omitted,  are 


1 30M) 

a(z,x) 


- 0(1(1)  = D2(M) 


(14) 


2 1 3 r.  2 /a  *-i  2D^(t)0(t) 

72  77— + -rr: T. 


q c(z,x) 


qc(1  - xc) 


26M 


1 + 


<?(1  - x*). 


= (1  - x2)—  + £2D^(ty) 

8x 


G(M) 


05) 


1 6(V,P)  p 

U _ = 'v’  P* 

q c(z,x) 

Using  the  values 

a = 6371  km 
h 


= 4 km 


g 

CL 


3 -1 

10  cm  sec 

7.29x10“^  sec”1 


p = 1 ,0  gm  cm 


-3 


= 2 dyne  cm 


-2 


M 

Si  = 

K = 
K_  = 


-7  2 -1 

10  gm  cm  sec 

2x10^  cm2  sec-1 

oc,  2 -1 

25  cm  sec  , 


(16) 


*•  t 


the  characteristic  dimensional  and  nondimensional  constants  are 
found  to  have  the  values 

= 1 .50x1 06  m'  sec’^ 

U = 6.89  cm  sec”"* 

Y = 1.60X10"3 

a = 6.28x1  o"*^ 

6 = 7 .i;2x10“^ 

£ = 8 ,5?x1 0-^ 

= 5.93xio"2 

X = 5.81 

K/K  = 800. 

r 

The  appropriate  values  of  K,  and  Kp  are  uncertain.  The  valuer 
chosen  above  are  based  on  estimates  given  in  Appendix  II. 

On  the  basis  of  the  study  of  relative  magnitudes  of  the 
terms  in  the  equations,  the  following  set  of  simplified  equations 
is  chosen  for  further  analysis: 

1 3(t,M)  P 

5 __ G(ljj)  = D (M) 

q d(z,x) 


G(M)  = (1  - x2) — -s-£_2D^(f)  ( 

3x 


i a(t,P)  p 

U = ?2p. 

q 3 ( z,x) 


This  choice  of  equations  involves  the  assumption  that  the  inequal- 
ities 


&<£2  \ — — DM2(t)]  « 1 

q 0(z,x) 


»#-  y 


- 15  - 


26M 


« 1 


hold  In  the  ocean  region  under  consideration. 


6.  Perturbation  equations 


As  the  system  of  equations  chosen  in  the  preceeding  section 
is  still  too  complex  to  treat  analytically,  the  equations  are  re- 
duced to  a simpler  form  by  the  standard  method  of  small  perturba- 
tions. Considering  a as  an  independent  parameter,  the  equations 
corresponding  to  the  limit  a -*■  0 are  determined.  These  equations, 
termed  the  zoro-ordor  equations,  are  independent  of  a and  have  the 
form 


2 » 
3?~ 


8M, 


- x 


as 


?,9p. 


- (1  - X ) 

ax 


(20) 


32m  at 

o . _ ro 

— 1 5 — + X ■ 1 — ■ 

az^  as 


= 6 


a(z,x) 


(21) 


92P, 


az 


2s  - P 


with  the  boundary  conditions 

to  = c)t0/c5z 


80  /oZ 
o 


and 


a(*0,P0) 

3 (z,x) 


= 0 


= 0 


M = 0 L at  z = 0, 


t = 0 


(22) 


a2t0/azr- 


(\4£)(i  - x2)^ 


m 


at  z 
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3M  /dz 

O’ 

Sp0/3z 


AT 

Q, 


at  z = 1 , 


For  a ^ 0f  correction  terms,  i.e*,  the  perturbation  terms, 
which  are  functions  of  a must  be  added  to  the  solutions  of  the 
zero-order  equations.  These  correction  terms  may  be  developed 
into  power  series  in  a.  As  a is  very  small  in  the  present  model, 
the  correction  terms  are  small  and  will  not  change  the  qualitative 
features  of  the  zero-order  solutions.  The  correction  terms  for 
a / 0 are  not  considered  in  the  subsequent  analysis. 

In  view  of  the  amount  of  simplification  required  to  arrive 
at  a system  of  equations  which  can  be  studied  analytically,  it  is 
of  considerable  Interest  to  know  which  processes  are  represented 
by  the  zero-order  equations  (20),  (21),  and  (22).  It  can  be  seen 
from  equation  (20)  that  the  zonal  flow  is  geos trophic  wherever  the 
term  €^8^f0/Sz^  is  negligible.  As  this  term  is  not  negligible 
near  the  upper  and  lower  boundaries  of  the  ocean,  the  zonal  flow 
hero  will  not  be  geos trophic.  The  process  of  eddy  diffusion  is 
effective  in  the  vertical  direction  only  because  the  vertical 
gradients  of  density  and  angular  momentum  are  much  larger  than 
the  horizontal  gradients.  The  zonal  component  of  wind  stress  is 
transmitted  almost  unchanged  to  the  ocean  bottom  because  the 
transport  of  zonal  angular  momentum  out  of  the  ocean  region  by 
the  meridional  currents  is  negligible.  The  nonlinear  term  repre- 
senting the  transport  of  relative  zonal  angular  momentum  by  the 
meridional  currents  is  included  in  (21 ) so  that  its  effect  on  the 
torque  acting  on  the  ocean  bottom  may  bo  determined.  The  transport 
of  mass  by  the  meridional  currents  is  significant  because  the  para- 
meter jj.  is  not  negligibly  small. 

In  the  subsequent  equations  the  subscripts  denoting  the 
zero-order  are  omitted  because  only  the  zero-order  equations  are 
considered. 


7.  Currents  driven  by  surface  wind  stress  in  a homogeneous  ocean 


The  simplest  flow  described  by  the  present  model  occurs 


when  the  ocean  has  a homogeneous  density  structure,  ise., 
mass  flux  across  the  ocean  surface  is  zero. 


v/hen  the 


I 


are 


The  equations  governing  the  motion  in  the  homogeneous  ocean 


2 3** 


8M 

x — = 0 

Sz 


(23 


at 

+ x — 


Bz 


a(t,M) 

6 . 

3(z,x) 


(2k 


These  equations  are  of  the  boundary- layer  type  because  the 
highest-ord9r  derivative  is  multiplied  by  a small  parameter,  A 
first  approximation  to  the  solution  of  (23)  and  (2ij.)  can  be  obtained 
by  neglecting  the  nonlinear  term.  The  linearized  equations  are 

e — r + = AxT  (25 

az4 


3M 

— + xt  = XT.  . (26 

3z 

The  boundary- layer  solution  of  (25)  and  (26)  in  the  upper 
boundary  layer  is 

t = (XT./x)£i  - exp L~k ( 1 - z)]  cos  k(1  - z^ 

M = AT/kjll  + ^exp£-  k(1  - z)]£  cos  ^(1  - z ) - sin  k(1  - 

and  in  the  lower  boundary  layer  is 

t = (XT/x)  - exp(-  kz)(  cos  kz  + sin  kz)[] 

M = (XT/k)[l  - exp(-  kz)  cos  kz}, 


where 
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In  the  interior  of  the  ocean,  the  solution  is 

t = XT/x 
M = >vT/k . 


As  long  as  the  upper  and  lower  boundary  layers  are  separated,  the 
boundary-layer  solution  may  be  expressed  compactly  as 

t = (XT/x)(1  - f)  (27) 

M = (XTA)O  - P),  (28) 

where 

f = expQ-  k(1  - z )3  cos  k ( 1 - z)  + exp(-  kz)(cos  kz  + sin  kz)  (29) 

3P/3Z  = - kf,  F(0,k)  = 1, 

The  current  system  given  by  (27)  and  (28)  is  the  same  as 
the  current  system,  set  up  by  a steady  wind  stress  as  described  by 
Ekman  ( 1 905 ) • 

By  substituting  (27)  and  (28)  into  (2 1+) , the  bottom 
torque  Tg,  including  the  contribution  of  the  nonlinear  term,  may 
be  computed.  The  bottom  torque,  to  the  first-order  approximation 
in  6,  is  given  by 


>6  d /T 


TB  = T + 


8k  dx\x 


and,  depending  on  the  distribution  of  T,  can  be  greater  or  less 

than  the  surface  torque.  As  6x/8k  is  very  small  in  the  present 
model,  the  bottom  torque  cannot  differ  appreciably  from  the  sur- 
face torque,  and  no  qualitative  changes  occur  if  the  nonlinear 
terms  are  neglected  altogether. 

The  current  system  set  up  by  a meridional  component  of 
wind  stress  is  confinod  to  the  upper  layers  of  the  ocean.  For 
simplicity,  'c  is  assumed  to  be  zero  in  the  subsequent  analysis. 


8.  Currents  driven  by  a mass  flux  across  the  ocean  surface 


If  the  surface  wind  stress  is  zero  everywhere,  the 
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circulation  set  up  by  a mass  flux  across  the  ocean  surface  must 
satisfy  the  equations 


^ o. 

, ov9P 

1 » : + x'ir  = 
az4 

- (1  - x-)— 
3x 

a2p 

. \ 

(d'ty  3p  3$  3p 

3z 

^\az  3x  ax  az 

3M 

— + xf  = 
3z 

0. 

) 


(31 ) 


(32) 


No  exact,  nontrivial  solutions  of  these  equations  have 
been  obtained  because  of  the  nonlinearity  of  (31).  The  formal 
boundary-layer  equations  are  too  complicated  to  solve.  I fc  is 
possible,  however,  to  obtain  approximate  solutions  which  are 
valid  in  the  interior-  of  the  ocean.  In  order  to  relate  the 
interior  density  distribution  to  the  surface  mass  flux,  the 
boundary- layer  equations  are  simplified  by  taking  an  average 
value  of  dl|r/3z  in  the  boundary  layer  and  by  neglecting  smaller 
terms  in  the  equations , 

The  boundary- layer  solution  of  (30)  is  given  with 
sufficient  accuracy  by 

1 - x2  3p 

* 5 (1  - t), 

x 3x 

where  f is  given  by  (29)0  This  boundary-layer  solution  for  $ is 
obtained  by  neglecting  the  vortical  variations  of  3p/dx  in  the 
boundary  layers.  Using  the  value  37T/J|.k  as  a measure  of  the 
thickness  of  the  upper  boundary  layer,  the  average  value  of 
3$/3z  in  the  upper  boundary  layer  is  estimated  to  bo 

a?  4k  i - *?  ap 

3z  37T  -J?  ax* 
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In  deriving  this  estimate  of  3^/3 z,  the  vertical  variations  of 
3p/3x  in  the  boundary  layer  are  neglected. 

The  simplified  boundary- layer  equations  for  the  density 
may  be  written 


(33) 


The  density  gradient  in  the  boundary  layer  is  estimated 
by  integrating  (33)  with  respect  to  z again  neglecting  the  ver- 
tical variations  of  3p/3x;  The  gradient  is 


8p 

3z 


Vk  1 

Q - — - 

37 r 


(1 


- z ) . 


Below  the  boundary  layer,  i.e.,  at  z - 1 - 3fl/4k,  the  density 
gradient  is 


ap 

3z 


and  the  density  is 


P 


P 


surface 


3 ir/u  1 - x2 1 

3p\ 

8k  x*”  1 

W 

31T 
~ Q. 
4k 


(34) 


(35) 


Obviously,  these  estimates  of  p and  3p/3z  below  the 
boundary  layer  are  rough  estimates  only.  In  view  of  the  uncer- 
tainty/ in  choosing  appropriate  values  of fx  and  k,  the  increased 
complexity  of  making  more  exact  estimates  of  p and  3p/3z  is 
probably  not  justified.  A more  detailed  analysis  of  the  boundary- 
layer  equations  shows  that  the  estimates  given  above  are  quali- 
tatively correct. 

The  relationships  given  by  (34)  and  (35)  are  used  to 
extend  the  interior  density  and  density-gradient  distributions 
continuously  to  the  surface.  Similar  relationships  are  used  in 
the  lower  boundary  layer. 

The  equation  governing  the  interior  density  distribution 
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is  derived  by  neglecting  the  term  €2a^$/3z^  in  (30)  and 
eliminating  the  stream  function  from  (31 )•  The  equation  is 


a2p 

pp 

-/ 

i - x2  ap' 

\ 1 - x2  3p  S2p  “ 

Bz^  1 

1*  — 
|_3z 

3xV 

x2  3x1 

r x*”  3x  ozd^ 

• 

As  general  solutions 

of 

(36)  have 

not  been  found,  the  mass  flux 

(36) 


the  surface  cannot  be  specified  arbitrarily  but  must  be  chosen  to 
be  consistent  with  the  interior  density  distribution  given  by 
special  solutions  of  (36). 

A new  coordinate  variable  y is  introduced  so  that 


i - x a 
x ax 


6y 


that  is. 


y = i log, 


■fci)- 


x. 


Equation  (36),  in  terms  of  y,  becomes 

-2_  _ .2. 


g(y) 


5,2 

6 P 


3z 


{dp  a2p  ap  a2p  \ 

? 1 

k6z  dy'  dy  3z3y / 


(37) 


where 


1 “ fx(y)]' 
8(y)  = - ■ TT2 


[x(y)] 

A solution  of  the  form 

p = o(z)^(y) 

is  substituted  into  (37)  and  the  resulting  equation  is  written  as 

d2o  /do  r &2x  /d*\2i  / 

sSz  U \Am  -r- 


where *V  is  a constant.  There  is  no  loss  in  generality  in  assuming 


22 


"V  =1,  so  that  the  equation  for  o becomes 


d2o 


dz 


ao 

= - 2ua  — . 

' dz 


(39) 


A solution  of  (39)  is 

a — 


cQ  cotan  yUcQ(z  - zq). 


where  zq,  cq  are  arbitrary  constants. 


The  equation  for  X is 

f&x\2 


a2* 


— 1 - X — *y  = 2g^. 

.ay/  dy^ 

No  exact,  nontrivial  solutions  of  (40)  have  been  found.  It  is 
possible,  however,  to  obtain  an  approximate  solution,  which  is 
valid  over  a limited  range  of  y,  of  the  form 


(40) 


where 


g, 


* = Me 


= exp (a  + a.y  + apy^). 


The  constants  aQ,  are  chosen  so  that  g/gQ  is  approximately 

equal  to  unity  in  the  region  under  consideration.  This  may  be 
done  with  considerable  accuracy  for  an  ocean  lying  poleward  of  45° 
because  g(y)  is  asymptotic  to  an  exponential  function  of  y for 
large  values  of  y.  The  equation  for  X0  is 

1 2 .2. 


/d*D\  ^ d 2 g 

— / -*o~2T  = 2Mo  + 2-*c 

\dy  / a:y 


53  2a2^o  + 2*o 


(41 ) 


The  first  integral  of  (41)  is 


[—)  = - ^a2*o  loge  *o  + 


(42) 
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where  is  an  arbitrary  constant.  Thus,  the  sol.ution  of  (41), 
in  inverse  form,  may  be  written 

fx° 

7 ~ ~ Va  (I  - loge|  + 01|2)'5d^, 

where  yQ  is  an  arbitrary  constant.  The  integral  must  be  evalu- 
ated numerically. 

An  approximate  solution  of  (37)  is,  therefore,  given  by 
p(z,y)  = - cogo^Q  cotan[^c0(z0  - z)]  . 

Only  z and  c of  the  undetermined  constants  are 
specified  by  the  boundary  conditions.  The  constants  aQ, 
depend  on  the  region  chosen  for  the  ocean,  and  the  constants 
y , c^  determine  a particular  horizontal  distribution  of  surface 
mass  flux.  The  details  involved  in  relating  the  boundary  condi- 
tions to  the  interior  density  distribution  are  given  in 
Appendix  III. 


9 . Currents  driven  by  both  surface  wind  stress  and  surface  mass 


If  the  circulation  is  influenced  by  wind  stress  and  sur- 
face mass  flux,  the  equations  governing  the  circulation  and  the 
density  distribution  have  the  form 

2 ^ ^ 'jl.  p 

€ — j-  + x"t  = XxT  - (1  - x ) — 

3z4  dx 


r[ 


3p 


3z  ox 


at  3p 

3x  3z 


3M 

— + xt  = XT. 

3z 


(43) 


(44) 


flux 


(45) 


(46) 


(47) 
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The  boundary-layer  solution  of  (45)  may  be  written 


- f). 


where  f is  given  by  (29)  . The  average  value  of  dty/dz  in  the 
upper  boundary  layer  is  approximately 


4k /XT  1 - x2  3p\ 

31U  x x dxj 

so  that  the  simplified  boundary-layer  equation  for  the  upper 
boundary  layer  is 


a2p  4/ik/xT  i - x2  sp\ap 

3z2  3 IT  \x  7 dxjdx 


(43) 


An  approximate  value  of  3p/3z  below  the  upper  boundary  layer  is 
found  by  integrating  (48)  with  respect  to  z neglecting  the  vari- 
ations of  ap/3x  in  the  boundary  layer.  The  value  is 


3p 

— fc? 

3z 


ap\ap 
a x/  ax 


(49) 


Elimination  of  the  stream  function  from  (46)  yields  the 
following  equation  for  the  interior  density  distribution t 


3p  a2p  ' 

3x  3z3x_ 


(50) 


As  the  term  Ad(T/x)/dx  is  an  order  of  magnitude  larger 
than  the  horizontal  density  gradient  in  the  interior  of  the  ocean, 
the  terms  involving  dp/dx  in  ($0)  are  neglected.  The  resulting 
simplified  equation  is  linear  and  a solution  may  be  chosen  of  the 
form 

p(z,x)  = p1  (x)^xpj^xj-- -j(zQ  - z)3  - d^,  (51  ) 


where  zo,  dQ  are  arbitrary  constants. 


r 
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The  interior  density  distribution,  given  by  (51),  is 
fitted  to  the  density  distribution  in  the  upper  boundary  layer 
by  requiring  that  the  equation 


Q (1  - dQ) 


XT 

— - (i  - d0>. 

x 


1 - x dp. 


dx 


dp-, 

3x 


(52) 


obtained  by  substituting  (51  ) into  ( lp9 ) and  setting  z = zq,  be 
satisfied  at  z = 1 - 371/1^. 

Equation  (52),  together  with  an  analogous  equation  for  the 
lower  boundary  layer,  is  solved  numerically  for  p^ . 

The  interior  density  distribution,  as  given  by  (51),  is 
stable  if  d(T/x)/dx  is  negative,  i.e.,  if  the  vertical  component 
of  velocity  induced  by  the  wind  torque  is  positive.  For  simplic- 
ity, the  wind  torque  is  assumed  to  be  of  the  form 


T = - x(b  + b,,x) 

so  that 


d(T/x)/dx  = - b1 , 

where  b,,  is  a positive  constant. 


10.  Numerical  examples 


The  particular  distributions  of  wind  stress  and  mass  flux 
chosen  in  the  previous  sections  are  not  intended  to  resemble  the 
actual  distribution  of  these  quantities  over  the  ocean  surface. 
The  wind  stress  and  mass  flux  distributions  are  chosen  to  simpli- 
fy the  analytical  treatment  and  the  numerical  computations  and  to 
give  a density  field  resembling  qualitatively  the  density  field 
associated  with  the  Antarctic  Circumpolar  Current  when  the  com- 
bined influence  of  wind  stress  and  mass  flux  is  considered. 

The  ocean  region  is  chosen  to  lie  between  70.1 °S  and 
i|5°S,  roughly  the  region  occupied  by  the  Antarctic  Circumpolar 
Current.  In  this  region,  an  example  of  each  of  the  following 
three  types  of  flow  is  considered  in  detail:  (A)  the  circulation 
In  a homogeneous  ocean  sot  up  by  a distribution  of  wind  stress 
over  the  ocean  surface,  (B)  the  circulation  set  up  by  a 
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distribution  of  surface  mass  flux  alone,  (C)  the  circulation  set 
up  by  the  combination  of  surface  wind  stress  and  surface  mass 
flux. 

Streamlines  of  the  meridional  circulation  Illustrating 
the  three  different  types  of  flow  are  shown  in  meridional  cross- 
sections  in  figures  (1),  (2),  and  (3).  The  streamlines  are  lab- 
eled in  units  of  At|7{  = 8.72  x 10^  m3  sec-"'  ). 

The  density  structure  in  figures  (2)  and  (3)  is  indicated 

by  isopycnals  (broken  lines)  which  are  labeled  in  units  of 

py(=  1*6  x 10  J gm  cm  J).  For  purposes  of  comparison,  an  example 

of  the  observed  density  structure  associated  with  the  Antarctic 

Circumpolar  Current  is  shown  in  figure  (I4.)  . Hero,  the  broken 

1 

lines  represent  lines  of  constant  ot. 

The  vertical  distributions  of  the  zonal  velocity  component 
for  the  three  types  of  circulation  are  compared  in  figure  (5)» 

The  comparison  is  made  at  45>°S  whore  the  zonal  velocity  component 
for  each  type  of  circulation  reaches  a maximum.  The  horizontal 
distributions  of  surface  wind  stress  and  surface  mass  flux  in  the 
computations  are  shown  in  figure  (6). 

The  total  volume  transport  through  the  meridional  cross- 
section  for  the  three  types  of  flow  is  found  to  be  as  follows: 


(A) 

112  x 

106 

-1 

sec  , 

for 

wind 

s tress 

alone 

(B) 

- 7 55  x 
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sec  , 
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mass 

flux  alone 

(C) 

433  x 

106 

m3 

-1 

sec  , 

for 

wind 

s tress 

and  mass  flux. 

1 1 . Discuss  ion 


The  present  model  consists  of  an  ocean  zone  that  is 
uniform  and  uninterrupted  in  the  longitudinal  direction,  hence, 
a toroid.  The  model  ocean  is  of  constant  depth,  is  bounded  below 
by  a level  bottom  and  above  by  a free  surface,  and  is  located 
poleward  of  45°S . The  model  is  open  toward  the  equator  and  the 


1 3 

The  quantity  is  equal  to  10  (p  - 1),  where  p represents  the 

density  of  sea  water  brought  isothermally  to  atmospheric  pressure. 
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pole.  It  Is  assumed  that  no  flow  occurs  across  the  poleward  side. 

The  circulation  is  driven  by  wind  stress  and  by  the  flux 
of  heat  and  water  across  the  ocean  surface.  Exchanges  of  heat  and 
water  across  the  surface  give  rise  to  changes  in  the  density 
through  thermal  expansion  and  the  concentration  or  dilution  of  the 
salts  dissolved  in  sea  water.  The  presence  of  density  differences 
in  the  ocean  influences  strongly  the  types  of  circulation  which 
can  exist.  The  thermodynamic  processes  which  influence  the  circu- 
lation are  taken  into  account  in  a simplified  manner.  The  volume 
changes  associated  with  the  density  changes  are  neglected.  The 
actual  fluxes  of  heat  and  water  are  replaced  in  the  model  with  an 
equivalent  net  flux  of  mass*  In  the  interior  of  the  ocean,  mass 
is  distributed  by  convection  and  diffusion. 

Simple  distributions  of  surface  wind  stress  and  surface 
mass  flux  are  chosen  and  the  steady  three-dimensional  circulation 
which  can  be  induced  is  illustrated  by  the  following  three 
examples:  (A)  the  circulation  set  up  by  a west  wind  alone,  (B)  the 
circulation  set  up  by  surface  mass  flux  alone,  (C)  the  circulation 
set  up  by  the  combination  of  wind  stress  and  surface  mass  flux. 

The  meridional  circulation  of  type  A,  resulting  from  an 
eastward  stress  on  the  surface,  is  illustrated  in  figure  1 , The 
ocean  is  homogeneous  in  density,  and  the  steady-state  meridional 
currents  are  confined  to  relatively  shallow  layers  near  the  upper 
and  lower  boundaries  of  the  ocean.  These  meridional  currents 
transport  a negligible  amount  of  angular  momentum  and  no  net  mass 
into  the  model,  so  that  a relatively  simple  meridional  circulation 
can  satisfy  the  conditions  necessary  for  a steady  state  to  exist. 
The  condition  of  zero  flow  across  the  poleward  limit  of  the  ocean 
region  is  satisfied  by  choosing  the  wind  stress  to  be  zero  at  this 
limit.  The  zonal  current  is  uniform  with  depth  except  in  the 
boundary  layers.  For  the  wind-stress  distribution  shown  in  figure 
6 the  eastward  volume  transport  of  the  zonal  current  is 
112  x 106  m3  sec"1  . 

If  there  is  a positive  upward  flux  of  mass  through  the 
ocean  surface,  nonhomogeneous  density  die tributiom can  exist  in 
the  ocean,  as  illustrated  in  figures  2 and  3*  The  mass  lost 


28  - 


through  the  ocean  surface  must  be  replaced  by  meridional  currents. 
In  the  circulations  of  type  B and  C,  the  mass  lost  through  the 
surface  is  replaced  by  a relatively  strong,  shallow,  poleward 
current  flowing  along  the  bottom  of  the  ocean.  In  type  B,  the 
density  decreases  along  streamlines  only  in  the  poleward  current, 
except  in  the  surface  layer,  while  in  type  C,  the  density 
decreases  along  streamlines  throughout  the  entire  region.  In  the 
absence  of  wind  stress,  the  mass  flux  induces  a strong  westward 
zonal  transport  totaling  755  * 10  nr  sec-  . The  presence  of 
eastward  wind  stress  increases  the  poleward  current  and,  conse- 
quently, the  mass  transport  in  the  bottom  layer  to  such  an  extent 
that  the  horizontal  gradient  of  density  throughout  most  of  the 
ocean  is  opposite  in  sign  to  the  gradient  induced  by  mass  flux 
alone.  The  total  eastward  transport  is  thus  larger  than  the 
transport  due  to  the  wind  stress  alone,  amounting  to  a total  of 
433  x 106  m3  sec"1. 

The  present  model  makes  possible  a more  quantitative 
estimate  of  the  Influence  of  thermodynamic  processes  on  oceanic 
circulation.  The  examples  of  the  circulation  discussed  here  are 
of  a very  special  nature  and  have  been  obtained  under  restrictive 
assumptions  so  that  quantitative  application  of  results  based  on 
these  examples  to  the  actual  ocean  is  not  justified.  The  examples 
are  given  to  illustrate  some  general  features  of  the  present  model. 

A characteristic  feature  of  the  present  model  Is  that 
whenever  the  eddy  flux  of  mass  or  of  angular  momentum  is  nondiver- 
gent,  the  meridional  component  of  velocity  is  negligibly  small. 
Thus,  if  the  present  model  can  be  validly  applied  to  the  actual 
ocean,  the  conclusion  Is  reached  that  any  steady  meridional  current 
must  be  accompanied  by  strong  diffusion  or  mixing  processes  which 
change  the  density  and  the  angular  momentum  of  the  water  along  a 
particular  streamline.  The  strongest  meridional  currents  develop, 
according  to  the  present  model,  near  the  boundaries  of  the  ocean 
where  the  eddy  stress,  i.e.,  the  momentum  flux,  can  attain  a 
relatively  high  value. 

The  total  zonal  transport  of  water  is  greatly  influenced 
by  the  presence  of  surface  mass  flux,  as  can  be  seen  by  comparing 
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the  total  zonal  transports  for  examples  A and  C.  Thus,  it  would 
be  Incorrect  to  explain  the  zonal  transport  of  the  Antarctic 
Circumpolar  Current,  which  is  computed  from  the  observed  density 
distribution,  in  terms  of  the  surface  wind  stress  alone.  Accord- 
ing to  the  present  model,  a considerable  range  of  values  of  the 
zonal  transport  is  obtained  for  the  same  wind  stress  by  varying 
the  distribution  of  surface  mass  flux. 

It  has  not  been  found  possible  to  justify  the  assumption 
that  processes  occurring  outside  of  the  ocean  region  under  consid- 
eration have  little  or  no  influence  on  the  circulation  in  the 
region  considered. 

The  solutions  obtained  are  not  valid  outside  of  the  region. 
Furthermore,  no  stable  density  distributions  have  been  found  for 
the  regions  where  there  is  a positive  downward  flux  of  mass  through 
the  ocean  surface. 

The  meridional  circulation  shown  in  the  diagram  differs 
considerably  from  the  circulation  believed  to  exist  in  the 
Antarctic  Circumpolar  Current.  The  difference  cannot  be  attributed 
to  oversimplified  boundary  conditions  used  in  computing  the  circu- 
lation, but  must  be  attributed  to  the  failure  of  the  model  to 
simulate  the  processes  occurring  In  the  actual  ocean. 
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Appendix  I,  The  mass -transfer  equation 

The  density  of  sea  water  Is  a function  of  essentially 
three  variables;  temperature  ©,  salinity  S,  and  pressure  p.  As 
the  dependence  of  density  on  pressure  is  not  dynamically  signif- 
icant, it  is  sufficient  for  the  purposes  of  the  present  study  to 
consider  the  density  as  a function  of  temperature  and  .salinity 
only  , 

The  heat  content  H and  the  salinity  S are  conservative 
quantities  in  the  ocean,  and  it  is  assumed  that  the  equations 
governing  their  distribution  are  of  the  form 

DH/Dt  = K9V72H 
BS/Dt  = Ks<72S, 

where  D/Dt  represents  the  total  time  derivative,  KQ  and  Kg  the 
constant,  isotropic  kinematic  eddy  coefficients  of  heat  conduc- 
tivity and  salt  diffusivity  respectively.  The  compressibility  of 
the  water  is  neglected  in  choosing  the  equations  to  be  of  this 
form. 

The  heat  content  H,  relative  to  an  arbitrary  reference 
point,  can  be  expressed  in  terms  of  temperature  by  the  formula 

h = cpo  - eo>, 

where  c^,  the  specific  heat,  is  assumed  to  be  constant. 

If  the  density  is  a linear  function  of  the  temperature 
and  salinity  of  the  form 

p = p[ 1 - a0(0  - 90)  + as(S  - SQ)], 

where  a0  and  a.g  are  constants,  and  if  the  eddy  coefficients 
Kq  and  Kg  are  equal,  the  two  conservation  equations  for  heat  and 
salt  can  be  combined  into  a single  equation  of  the  form 

Dp/Dt  = Kp<72p, 

which  represents  the  conservation  of  mass  under  convective  and 
diffusive  processes. 
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The  actual  relationship  between  density,  temperature, 
and  salinity  may  be  approximated  by 

P = P0D  - ae<e  “ ©0>  + aS(S  ' So0» 

where  <Xq  and  as  are  functions  of  both  temperature  and  salinity, 

r.na  where  p , 9 , and  S are  constant  reference  values.  The 
ro ’o’  o 

coefficients  ctg,  can  be  calculated  from  tables  . Using  the 
auxiliary  condition 

Gq(9,S)  - q.q(Q f& ag(9,S)  - ag(9^,S) 

S - 9-9 

o o 

and  the  reference  values  9 = 10  C,  S = 35%,,  and  p = 1 .027, 

0 ° ° _i. 

the  coefficient  <Xq  was  found  to  vary  between  0.96  x 10  4 C 
and  2,06  x 10-^  for  temperature  between  -2°C  and  18°C  and 
salinities  between  32$oand  37%o.  Over  the  same  range  of  temp- 
erature and  salinity  the  coefficient  a.g  was  found  to  vary 
between  0.750  and  0.776, 

The  values  of  <Xq  and  ag  used  in  the  present  report  are 
<Zq  • = 1 .00  x 10"^  C"1 

ctg  ~ 0.76* 

The  flux  of  mass  across  the  ocean  surface  may  be  expressed 
in  terms  of  the  flux  of  heat  and  the  flux  of  water  across  the 
su^ft'/vCc  The  flux  of  water  is  equal  to  the  difference  between 
the  rates  of  evaporation  E and  the  precipitation  P.  If  E - P is 
negative,  i.e.,  water  is  added  to  the  ocean,  the  salt  content  is 
decreased  in  the  surface  layers.  This  decrease  in  salinity  is 
interpreted  as  being  due  to  a virtual  flux  of  salt  Qg  across  the 
ocean  surface  equal  to  pS(E  - P),  wher-8  S is  the  mean  surface 
salinity. 

The  boundary  condition  for  the  density  at  the  ocean 
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surface  may,  therefore,  be  written 


where  Q,  Q,q,  and  Q,g  are  the  fluxes  of  mass,  heat,  and  salt 
respectively. 

It  is  difficult  to  estimate  with  any  degree  of  accuracy 
the  values  of  and  Qg,  Generally  speaking,  the  ocean  loses  heat 
to  the  atmosphere  in  the  region  of  the  Antarctic  Circumpolar 
Current  but,  because  the  precipitation  exceeds  evaporation,  the 
density  at  the  ocean  surface  remains  lower  than  the  density  at 
greater  depths.  Assuming  that  the  maximum  net  precipitation  rate 
is  of  the  order  of  130  cm  year  and  the  minimum  net  heat  loss  is 
of  the  order  of  .10  cal  cm”  day"  , and  using  the  average  surface 
salinity  of  3I4 the  values  of  Qq  and  Qg  are  found  to  be 

Qq  = - 1 ,2  x 10”^  cal  cm”^  sec""* 

“7  • 2 

Qg  = - 1 ,4  x 10  ' gm  cm"  sec"  . 

Thus,  an  estimate  of  the  maximum  mass  flux  according  to  the 
definition  given  above,  is 

QM  = - (ae/cp)Q0  + ctgQg  *=■  10“'  gm  cm"  sec"  , 

This  value  of  is  assumed  to  be  representative  of  the  actual 
maximum  mass  flux  across  the  ocean  surface. 


Appendix  II . The  kinematic  eddy  coefficients 

The  boundary  conditions  in  nondimens ional  form  for  the 
zonal  angular  momentum  and  the  density  at  the  ocean  surface  are 
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where  and  are  maximum  values  of  the  wind  stress  and  mass 
flux  respectively. 

An  estimate  of  Q..,  given  in  appendix  I,  is 

-7  -2  -1  M -2 

10  gm  cm  sec  . The  value  of  2 dyne  cm  is  chosen  for 

a3  being  appropriate  for  the  ocean  region  containing  the 

Antarctic  Circumpolar  Current. 

In  order  to  have  density  differences  in  the  present 
model  comparable  with  the  observed  density  differences  in  the 
region  of  the  Antarctic  Circumpolar  Current,  the  characteristic 
fractional  range  of  density  variation  y must  have  a value 
comparable  with  the  value  obtained  in  the  actual  ocean.  In  the 
present  model,  y is  given  b?/-  hQj^/pKp.  Thus,  for  = 25  cm  sec 

the  value  y is  found  to  be  1 .60  x 10’*^.  This  value  of  y is  of 
the  same  magnitude  as  the  fractional  range  of  density  variation 
observed  in  the  ocean,  A higher-  value  of  K , e.g.,  = 

2-1  P P 

100  cm  sec  , would  give  almost  negligible  density  difference 
according  to  the  present  model. 

The  eddy  viscosity  K is  chosen  to  have  the  value 
2 x 10^  cm^  sec"^  , A smaller  value  of  K would  give  excessive 
velocities  throughout  the  ocean.  A larger  value  of  K would  give 
smaller  velocities  but,  would  also  give  the  result  that  the 
circulation  is  determined  primarily  b^  the  mass  flux  across  the 
ocean  surface.  This  does  not  appear  to  be  the  case  in  the 
actual  ocean. 

The  necessity  of  assigning  a large  value  to  the  ratio 
K/Kp  in  order  to  obtain  reasonable  results  from  the  present 
model  reflects  the  lack  of  a more  suitable  mechanism  to  counter- 
act the  wind  torque  applied  at  the  ocean  surface. 

Appendix  III.  Boundary  conditions  (B) 


The  approximate  equations  relating  the  interior  density 
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gradients  to  the  density  gradients  in  the  boundary  layers  are 


~ b Q - - at  z = 1 - 3‘fp'l+k, 

3z  gldy/ 


ap  n ap 

5z  gldy 


at  z = 7T/k. 


These  equations  determine  the  magnitude  of  the  arbitrary 
constants  in  the  solution  for  the  interior  density  distribution. 
The  derivatives  3p/3z,  3p/3y  iiiay  be  written 


- co[cotan^c0(zo  - z)  + tan^cQ(zo  - z)]p 


— — P,  gn  = exp  (a  + a,y  + a2yd), 

ay  sQZo  dy 


so  that  at  z = fl/k, 


c0Ccotan/^co -zo  “ ^A)  + tanttc o(zQ  ~TT/k)]p  = - --^g p2. 

ss0 \ Gy  f 


Substitution  for  p yields  the  equation 


tan  " uc  (z  - 77/k)  = - 1. 

' gSo^0Vdy  / 


This  equation  is  satisfied  approximately  by  choosing  the 
arbitrary  constant  c.,  of  (42)  so  that  the  variation  of 

O * 


P 

GKgX,  )/dy|  /ggo>Q  is  minimized,  A relationship  between  zq 


and  c is  then  obtained  using  the  average  value  of 


[<l(g*0)/<lyj  /ggQ*0  in  the  region  under  consideration. 

The  mass  flux  Q,  across  the  ocean  surface  is  determined 
from  the  equation  connecting  the  interior  density  gradients  to 
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the  density  gradients  in  the  upper  boundary  layer,  i.e«. 


dp  u(  ap\2-  3ir 

Q = — + Cl— \ , z = 1 - — . 

bz  g V3y/  4k 

Substitution  for  p yields 


Q = •"  t 


^o\oco 


‘-'l 


co  tan  jUco 


(zo  “ 1 


3^410 


2 


( L cotan  ^c0(z0  - TT/k) 


The  recuiremert  that  Q has  a unit  negative  extremum  value  yields 

a second  relation  between  z and  c so  that  both  constants  can 

o o 

be  determined. 
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